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ON THE LOGARITHMIC -SINGULARITY CORRECTION IN 
THE KERNEL FUNCTION METHOD OF SUBSONIC 
LIFTING-SURFACE THEORY 

C. Edward Lan* and John E. Lamar 
Langley Research Center 

SUMMARY 

A new logarithmic-singularity correction factor is derived for use in kernel func- 
tion methods associated with Multhopp's subsonic lifting-surface theory. Because of the 
form of the factor, a relation has been formulated between the numbers of chordwise and 
spanwise control points needed for good accuracy. This formulation is developed and 
discussed. Numerical results are given to show the improvement of the computation 
with the new correction factor. 


INTRODUCTION 

There are various kinds of singularities which arise in the kernel function methods 
of the lifting- surface theory. Since Multhopp published his method in 1950 (ref. 1), much 
effort has been expended in obtaining an accurate accounting of these singularities. The 
singularities are associated with: (1) the mathematical modeling of the lift potential of 
a wing surface (i.e., doublets) and the resulting downwash produced by them as calculated 
in the analytical or numerical spanwise integration over the lifting surface, (2) leading- 
edge pressures, and (3) changes in wing sweep at the apex or along the span (refs. 2 and 3). 

Singularities associated with wii^ sweep are not given special treatment (i.e., 
special pressure modes), but instead Multhopp's procedure for rounding the leading and 
trailing edges in the vicinity of the affected regions is used (ref. 1). Leading-edge pres- 
sure singularities are expected because of the modes assumed, but they add to the first 
singularity problem listed and are treated in the context of this problem. The first singu- 
larity is the main subject considered herein. 

The spanwise integration over the lifting surface generally has associated with it a 
logarithmic -type singularity term at some spanwise position which in a purely numerical 
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treatment can lead to significant errors. To minimize these errors, the usual procedure 
is to separate the logarithmic term and evaluate it to provide a "correction" to the 
remainder of the integral, (Refs. 1 and 4 are examples.) The term which results is 
called the logarithmic-singularity correction (LSC). This approach was introduced by 
Multhopp in reference 1 (used later by others in refs. 2 to 9) along with a procedure for 
the LSC calculation. Improvements in the computational procedure have been formu- 
lated by Mangier and Spencer (ref. 9) and Zandbergen, Labrujere, and Wouters (ref. 2). 
However, the most commonly used correction has a serious deficiency which arises 
because the correction term consists only of the leading term of a correction series 
which diverges as the chordwise control points approach the leading and trailing edges. 
The exclusion of the high-order correction terms, logarithmic in nature, in this diver- 
gent series has been credited by Jordan (ref. 8) as the main cause of significant errors 
in some existing lifting-surface methods when the number of chordwise control points is 
increased. Jordan (ref. 8), Wagner (ref. 5), and Medan (ref. 3) have advocated the use 
of wing-edge control points for better accuracy in the solution. These wing-edge control 
points are known to be associated with algebraic, instead of logarithmic, singularity cor- 
rections. However, even if wing-edge control points are used, the logarithmic correction 
factors which are needed for interior control points are still inaccurate near the edges. 

Another feature of existing kernel function methods is that the relation between 
numbers of chordwise and spanwise control points used to obtain a reliable solution has 
been empirical, derived by applying the theory to a limited number of planforms (ref. 7). 
However, experience indicates that this empirical relation is not applicable to arbitrary 
combinations of configurations and Mach number. It is in fact highly dependent on plan- 
form and Mach number. 

In this report, a method of obtaining a convergent correction series is presented 
and the leading term of such a series is derived. Better convergence characteristics of 
a modified Multhopp' s method are illustrated when the new correction factor is used. In 
addition, a general relation governing numbers of chordwise and spanwise control points 
is derived and demonstrated to provide accurate solutions for arbitrary planar planforms. 

Some details of the development are given in appendixes A and B, and implementa- 
tion of the analysis is given in appendix C. 

SYMBOLS 


A 



-X) 


2 


+ Y 
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b ,b 
vv’ vn 




^D,i 


"Q! 


ACr 


c(i?),c(y) 


wing span, equal to 2 in this report 

Multhopp's quadrature weighting factors defined in equation (B21) 
ratio of near -field to far -field vortex drag 

overall wing lift-curve slope 

lifting pressure coefficient 
constant to be determined 

streamwise half -chord at rj and y, respectively, nondimensionalized 
by b/2 


E(k) 

E(<^,k) 

F(</),k) 

Gj(0) 

^jpv’^jvn 

K(k) 


average chord 

complete elliptic integral of second kind 
incomplete elliptic integral of second kind 
incomplete elliptic integral of first kind 
function of <p defined by equation (21) 
jth chordwise loading function 
jth chordwise integral defined by equation (6) 
see equation (B22) 

integrals defined by equations (15) and (16) 
complete elliptic integral of first kind 

modulus of elliptic integrals squared, ^ ~ ~ B)^ 



complementary modulus of elliptic integrals squared, 


1 - k2 - (A + B)2 - 1 
4AB 

see equation (B20) 

k' at left and right wing tip, respectively 
logarithmic component of Ij 
free-stream Mach number 

number of span stations where pressure modes are defined 
number of chordwise control points at each of m span stations 
jth spanwise loading fimction 

perturbation velocity in vertical direction nondimensionalized with respect 
to free-stream velocity 

2 civ) 

_ ^ - ^le^) 

2 c(y) 


^ - ^le(^) 
2 c(r?) 


= /3 


y - V 

2 civ) 


= ^ 


y - V 

2 c(y) 


streamwise and spanwise rectangular Cartesian coordinates nondimensiona- 
lized with respect to b/2 and associated with control-point locations 

aerod 3 mamic center in fractions of reference chord and referenced to leading 
edge of reference chord, positive aft 


/3 
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n(Q!2,k) 
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<0 

Subscripts: 

le 

n 

s 

te 

V 


= \jl - 

small positive number 

angle for locating control points along span, cos“^(r 7 ) 
sweep angle, deg 

Heuman^s lambda function, equation (B3) 

taper ratio, Tip chord index in equation (23) 

’ Root chord 

rectangular Cartesian coordinates, nondimensionalized with respect to b/2, 
in X- and y-directions, respectively, and associated with pressure doublet 
locations 

complete elliptic integral of third kind 

incomplete elliptic integral of third kind 

angular coordinate defined by X = - (1 - cos 0) 

2 

angular coordinate defined by Xi = - cos 

= sin"l(A - B) 


leading edge 

spanwise influencing station associated with particular 77 values 
chordwise control point associated with particular 0 values 
trailing edge 

spanwise control point associated with particular y values 
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Abbreviation: 


LSC logarithmic -singularity correction term 


The symbol 


denotes finite -part integration. 


MATHEMATICAL DEVELOPMENT 


In linearized steady subsonic flow theory for a thin wing, the lifting pressure ACp 
is related to the downwash distribution on the wing surface through the following singular 
integral equation: 


sill. 


^le(’7) (y - 


ACp(|,7?) 


1 + 


xr.i. 


[ ^(x-^)2+^2(y _^)2J 


d| dTj 


( 1 ) 


where the spanwise integration is defined by the finite-part concept. By introducing the 
new variables, 


X = 


2 c{v) 


Xi = 


2 c(tj) 


equation (1) becomes 




/ 


X - Xi 


1 + — 

^ y(X - Xj ^ + ^^y - 7/)/| c(7?jj 


c{ri) dX^ d?7 


In equation (3), the lifting pressure coefficient is usually approximated by 


N-1 


. 2 V 


j=0 


c(j?) 


( 2 ) 


(3) 


(4a) 
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where 


and 




COS j(p^ + cos (j + 
sin 


(4b) 


Xi = |(l-cos 


<t>i) 


It follows that equation (3) can be written as 


(4c) 


w(X,y) = i 


47T J -1 (y - n)^ 


N-l 

2 qj(7?) Ij(X,Y,r?) 
j=0 


cbj 


(5) 


where 


I.(X,Y,77) = I r 
1 ^ 'JQ 


2 pi cos + cos (j + l) 0 j^ 


sin (p^ 


X - X, 


^(X-Xj)2..Y^ 


dXi 


( 6 ) 


and 


Y = 


2 c(j7) 


Now, it has been shown (ref. 1) that Ij behaves like (y - 77 )^ln|y - 77 ] near ri = y. The 
factor (y - tj)^ is canceled by the same factor in the denominator of equation (5), so that 
a logarithmic singularity at rj = y appears. Accurate spanwise integration can be done 
only if this logarithmic singularity is carefully accounted for. 

Mangier and Spencer (ref. 9) removed the logarithmic component of Ij, called Lj, 
before performing the spanwise integration. This approach, also used here, leads to 
(eq. (B14)) 


yx,Y,77) = Ij*(X,Y,77) + Lj 

so that the term in the summation in equation (5) can be replaced with (eq. (B15)) 


(7) 
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( 8 ) 


q.(77) i^{x,Y,n) = qj{y)(ij - + [qjC’?) - qj(y^ij + qj(y) V 

Obviously, equation (8) involves Lj and therefore requires attention. As — y, 

L*(X,Y, 77 ) varies in X like L(X,Y,y) and thus 
J J 

Logarithmic ,■ f ^ 

Lj = component |lj(X,Y,i 7 ) - IjpJ,Y,y) = AW (9) 



r 

1 



-| 


(M 1 

II 

<1 

1 ^X cos + cos (j + 1)0J 

-1 + — 

X - Xj 
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->1 cos + cos (j + 1)0J^ 
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Xj - X 

dxA 


'x sin 


\l^ 

- Xj)2 + Y^ 



>w 

V ' 

1 

J 


It is observed that for small Y, the factors inside the brackets in this equation are small 
when Xj^ is not close to X and have significant magnitude only near Xj^ = X. For this 
reason, Multhopp (ref. 1) argued that for deriving the logarithmic correction term, it is 
reasonable to develop the function hj (eq. (4b)) into a Taylor series about the point 
Xj^ = X, or equivalently about (p-^ = <p. The result can be integrated exactly and then 
developed into a series for small Y to obtain the logarithmic components. Now, it is 
obvious that hj does not have a Taylor series expansion about the integration end 
points, Xj^ = 0 and Xj^ = 1. For example, by using equation (4c), the Taylor series 
expansion for hQ is obtained as 


ho(Xi) = ho(X) 


dh(j(X) 


dX. 


, d"“hn(X) « 

(X. - X) + i ^ (X. - X)^ + 

( 1 } 2 dX,2 V 1 ) 


_ ^^1/2 Xj - X (3 - 4X)(X^ - X)^ 

X^/2 2xV2(i _ x)l/2 ^ 8X^/2 (1 - X)3/2 


( 10 ) 


This series obviously does not converge near X = 0 and X = 1. This n onconvergence 
of the series expansion is due to the presence of sin 4>-^ = 2 JXj^^l - Xj^^ in the 
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denominator of hj. Therefore, a natural way to improve the series expansion is to 
expand only cos + cos (j + l)0j^ as follows: 

cos ](f)^ + cos (j + = cos + cos (j + 1)<^ + ^^^|cOS j(j)^ 


+ cos (j 




(Xi - X) + o(Xi - x)= 


4>i=(p 


= gi(0) + 52^*^) (^1 - X) + 0(Xj - X)' 


where 


( 11 ) 


gj(0) = cos j(f> + cos (j + 1)0 


( 12 ) 


and 


g (0) = -2 3 sin ]0 + (j + 1) sin (j + 1)0 
^ sin 0 


(13) 


Hence, for small Y, equation (9) can be written as 


Ai. 

J 7T Jq 


^1 gj(0) +g2(0)(Xj - X) 

X - Xi 

1 + ^ 

Jo ^Xi(l - Xj) 

V^(X - Xi)^ + 


dX^ - l](X,Y,y) 




1 + 


X - X, 


\lf^) 




dX, 


gg(0) ^1 -Xi 
^ ^0 M Xj 


X - X, 


1 + ^ — 


dXj - lUx,Y,y) 


(|y| « l) (14) 
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} 

I 


where equation (6) has been used and lj(X,Y,y) is the integral l^(X,Y,y) with the 
loading function expanded in accordance with equation (11). Obviously, the logarithmic 
correction term comes only from the integrals expressed in equation (14). Note that if 
higher order correction terms are desired, more terms should be retained in equa- 
tion (11), The two integrals of equation (14) can be integrated exactly to give 



= 7r(l - Ag(;^,k[) (16) 


Details of the integration are given in appendix A. 

In appendix B, the LSC is derived; a summary of the derivation is given here. The 
relation between Y and k' is shown to be (eq. (B2)) 

Y^ = 4k'^X^(l - X)^ 1^1 + k’^(3 - 8X(1 - X^J + o(k’®) (17) 

Equation (17) shows that small k' implies always small Y. However, small Y does 
not always imply small k’, depending on the value of X. It follows that it is more 
natural to develop and J 2 for small k’, rather than for small Y, to find the loga- 
rithmic terms. This choice represents a deviation from Multhopp's original method. In 
developing equations (15) and (16) for small k’, it is foimd that the logarithmic compo- 
nents are (eqs. (Bll) and (BIO)) 
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_sign (y - V) ^ (y - r?)2 ^ In k’ + o(k'3 in k') (18) 

2 c(y) sin cp «77 

-sign {y -<!)-§-:■ (y - ^ 1" >=’ + >" >=') 

c(y) sm (jf) dr] 

It follows from equations (14) to (16) and (18) to (19) that for small Y the dominating 
logarithmic component of Ij can be written as (eq. (B12)) 

L. = -sign (y - T?) I — GA(p) (20) 

] IT c(v) J 


Logarithmic 
component 
of J, 


Logarithmic 
component 
nf .T» 


where (eq. (B13)) 

G.(0) = cos 0 -I- j sin j(^ -t- (j H- 1) sin (j -t- 1)<^ (21) 

3 sin 0 

It is shown in appendix B that equation (20) can be reduced to the conventional form by 
expanding dk'/d ?7 and k' for small Y. 

From the preceding development the first term of equation (8) involves Lj and 
contributes the following result to equation (5) (eq. (Bl6)): 


-27T W I 


■')part 1 ' ■ I * ‘=2 ■ ’‘ 2 ) 


( 22 ) 


where k'^ and k^ are the values of k’ at the left and the right wing tip, respectively. 
The notation has been made compatible with Multhopp's notation (ref. 1). Reference 7 
gives a detailed explanation of the subscripts which have been introduced. For the second 
and third terms of equation (8), Multhopp’s method of integration can be applied by using 
the trigonometric interpolation formula of the type 


m-1 


. 

3 3 m + 1 


^ m 

X=1 


sin \9 


(23) 
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where 6 = and m is the number of spanwise integration stations over the 

n 2 m + 1 

whole wing. The integration of the second and third terms of equation (8) in equation (5) 
results in (eq. (Bl9)) 



where the prime on the summation sign means that the term with n = v is omitted from 
the sixmmation. Combining equations (22) and (24) reduces equation (5) to the following 
equivalent algebraic equation: 


-2rrw(,^s’yJ 



- kj + kg In kg 




sin 9^ dk'^ 

m + 1 dJ] 



-I ‘’mOjC’n) 



(25) 


where b and b are Multhopp's quadrature weighting factors and are defined in 
appendix B for convenience. The detailed expression for dik^jda] is also given in 
appendix B. Equation (25) represents the new system of equations for determining the 
values of the loading fimction qj. 


REGION OF VALIDITY FOR THE LSC 


From the definition, it is known that 


j^,2 ^ (A + B)^ - 1 _ (1 - X)^ + 2Y^ - 1 + 2AB 

4AB “ 


4^(1 - X)^ + 


+ Y‘ 


- X + +\/(l - X)^ + 

2^(1 -X)^ + Y^\/x^ +Y^ 


(26) 
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As indicated previously, the logarithmic-singularity correction term (LSC) must be 
removed before accurate spanwise integration across the second-order singularity can 
be performed. However, the LSC is only important if k' is small. In accordance with 
equation (17), small k’ implies small Y, say Y = 0(e), where e is a small positive 
number. To see how this requirement restricts the values of X, consider X = 0(e) 
also. Examination of equation (26) shows that 


k’2 = = 0(1) 

e 0(1) 

It follows that the LSC is not important in this case. In fact, references 5 and 8 show 
that the correction term is algebraic, O^jYQj^'^^j, at the leading edge rather than 
logarithmic. When combined with the second-order singularity in equation (5), a singu- 
larity, 0^jYQ|”^^^j, remains in the spanwise integration which must be performed with 

the finite-part concept. Since X = 0(e) implies that a control point is close to the wing 
leadir^ edge, there must be a transition region near the leading edge in which the loga- 
rithmic singularity changes to an algebraic one. In this transition region, any use of the 
LSC would be improper, as it would misrepresent the type of singularity. 

On the other hand, consider Y = 0(e), X = 0(1), and (1 - X) = 0(1). Then, by 
binomial expansion. 



0(e2)(i + 0(1) + 0(e2)) 
2|x||l -X|[l +0(e2)3 


where it is assumed that 0<X<1. Therefore, k’ is 0(e), and the LSC is proper in 
this case. In other words, to use the LSC properly, the chordwise control points should 
be located so that 

|x|»1y 1 Y = 0(e) (28a) 


13 



for control points near the leading edge and 

|(1-X)|»|y 1 Y = 0(e) (28b) 

for control points near the trailing edge. 

To analyze these conditions further, consider relation (28a). It is convenient to 
express both X and Y in terms of Yq. Using the first two terms of a Taylor series 
and the binomial expansion leads to 


2 c(? 7 ) 


1 X - 4jg(y) - ijg(y) iv - y) 



= Xf) + — 
0 ^ 


(1 - Xq) tan Ajg + Xq tan 


^ ^ "^le) [(^ ■ ^o) ^le 


+ Xq tan A^g 





(29) 


where 


l(.a„ V - tan A.,) = ^ 


<ly 




tan Ajg = 


dy 


tan A^g = 


dy 


(30) 


and 


Y = ~ 

0 2 c(y) 


(31) 
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Similarly, 




2 c(t 7) 2 c(y) 


1 + (y - v) 
c(y) 




Y ^ / \ 

= Yq + -^(tan - tan + ©(Yq^) 


( 32 ) 


p 

If terms of the order of Yq are neglected in equations (29) and (32), the condition for 
the validity of the LSC near the leading edge (eq. (28a)) becomes 


Xo + — 

0 (3 


Xq) tan Ajg + Xq tan A^^ 


» Y, 


(33) 


A critical test of relation (33) would be when Yq is negative, |y - t?| is as small as 
possible to ensure the importance of the LSC, such as in the tip region, and Xq is for 
the control point nearest the leading edge. In numerical calculation, | y - ^?| is small if 
the number of spanwise integration points, m, is reasonably large. Reasonably large m 
is necessary for accurate integration. Therefore, relation (33) provides a criterion to 
determine how close to the leading edge the control points are allowed to be located. 
Relation (33) shows that the condition involves the effects of aspect ratio and taper ratio 
(throi^h c(y)), Mach number, and the sweep angle. For example, relation (33) shows 
that the larger the sweep angle and subsonic Mach number, the more difficult it is to sat- 
isfy the relation. Most existing kernel function methods require empirical relations 
between the numbers of chordwise and spanwise control points for convergence. Such 
relations are usually applicable only to a limited set of planforms and a limited range of 
flow conditions. Relation (33) offers a theoretical relation to ensure the correct repre- 
sentation of logarithmic singularity and therefore to ensure better accuracy of the 
numerical integration. 

A similar criterion can also be set up for control points near the trailing edge. 
However, it is known that the behavior of Ij at the trailing ecfee in the spanwise integral 

is algebraic, being 0^| YQj^'^^j, compared with 0^|YQ|^/^j at the leading edge (ref. 5). 
Therefore, near the trailing edge the singularity in the spanwise integration of equa- 
tion (5) is at most 0^|YQ|“^'^^j which is integrable, even if the behavior in the inte- 
grals Ij has not been accounted for at all. The situation is quite different at the leading 
edge. Hence, a similar criterion for control points near the trailing edge is less impor- 
tant than that of relation (33) in the overall accuracy of computation. 
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CRITERION FOR SELECTION OF N,m SET 


In order to satisfy relation (33), reexpressing it in the following slightly different 
form is useful: 




Xq + — tan Ajg - — Xq (tan - tan A^^) 


» Y, 


( 34 ) 


Near the leading edge the critical test mentioned previously is for Yq becoming a small 
negative number. Thus, since A^g is in general larger than A^g, the smallest value of 
the left-hand side would occur if the third term were omitted because its contribution 
would be only a small positive number under these conditions. So, let 


Xo + tan Ai^ 


» 


(35) 


The left-hand side of relation (35) must be positive without the absolute value signs, 
since it represents X. (See the discussion following eq. (27).) Therefore, the following 
relation must be satisfied: 




tan A 


le 


(36) 


Cj^ Dependence 

In order to determine the magnitude of the terms required to satisfy the "much 
larger than" stipulations, it is expedient to reformulate relation (35) into a simple 
"larger than" test. This can be done by rewriting relation (35) as 


^0 ^le ^ ^1 


(37) 


with Cj being an as yet undetermined constant. In order for Xq to be as small as 
possible, it should be written for the first control point to give 



X - ^le^y) 

2 c(y) 


1 - cos 


2it 

2N + 1 
2 


(38) 
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where N is the number of control points in the chordwise direction. Similarly, for Yq 
to be as small as possible, it should be written for the integrating station nearest the wing 
tip and the control point at the next inboard location: 


_ j3(y -n) - 

0 " 2 



(m - 1)7T 
2(m + 1) 


sin 


(m - 3)tP /b\ 

2(m + li, \2/ 


(39) 


where m is the number of spanwise integration stations from tip to tip. Thus, the prob- 
lem of finding an N,m set reduces to determining a proper value for Cj^. The approach 
taken is to select different values of and then examine the consequences in terms of 
N and m which satisfy relation (37) and yield ratios of near-field to far-field vortex 
drag Cjj ^ near unity. According to Multhopp (ref. 1), the drag ratio information 

provides an indication as to the accuracy of the solution. This ratio, which ideally would 
have a value of unity, is very sensitive to variations in N and m and is obtained from 
the equations given in reference 7. The results of this study are presented in table I for 
three thin flat wings at M = 0. Wing 1 is an aspect -ratio -7 rectangular wing, wing 2 is 
an aspect-ratio-3.5 sheared rectangular wing with A = 45°, and wing 3 is an aspect- 
ratio-2 delta wing. 

In order to determine a satisfactory N,m set, initial values are needed. The 
initial values selected for all three wings were N = 8 and m = 23. The sets which 
satisfy relation (37) are listed in table I. The solutions for wings 1 and 2 indicate that 
Cj “ 10 would be reasonable. It is also clear from table I that (1) though N and m 
may satisfy relation (37), there is no guarantee that the Cj-j i value will be near 1; 

and (2) for wing 3, no satisfactory value of is determined. That wing 3 should pre- 

sent a convergence problem is perhaps not surprising, since the leading- and trailing- 
edge sweep angles of this configuration are markedly different, and thereby the term 
omitted in going from relation (34) to (35) is emphasized. (See the discussion which 
follows relation (34).) However, when relation (34) was used to determine the N,m set 
for wing 3 at = 10, no change in the set and consequently in Cp . occurred. 

(See table I.) The difficulty seems to be in the initial estimates of N and m. 

The slow convergence for wing 3 compared with that for the other wings is caused 
by a difference in Yq. For the other wings Yq is 0(e), whereas for wing 3, Yq is 
0(1). This change in the order of magnitude of Yq is caused by the diminishing chord 
at the tip of wing 3 and can be seen from an examination of equation (39). Note that from 
equation (38), Xq is 0(1) for all wings. From the preceding discussion, it is apparent 
that for wing 3 the left-hand side of relation (37) is the difference between two 0(1) terms 
thus, the relationship has less impact on a proper choice of N and m, or it becomes a 
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TABLE I.- EFFECT OF Cj ON N, m, AND ON ^ 

[m = o3 


/-I 

Wing 1 


Wing 2 


Wing 3 

^1 

N 

m 

^D,ii/^D,i 

N 

m 

^D,ii/^D,i 

N 

m 

^D,ii/^D 

2 

7 

49 

0.5105 

8 

47 

-0.3590 

8 

43 

-1.5514 

4 

6 

63 

-.1250 

7 

55 

.2571 

7 

47 

-.8792 

6 

5 

63 

.8997 

6 

57 

.7224 

7 

55 

-.5417 

8 

5 

73 

.9700 

6 

65 

.8687 

6 

51 

-.1010 

10 

4 

67 

1.0746 

5 

59 

1.0380 

6 

57 

^1055 

^10 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

12 

4 

75 

1.0883 

5 

67 

1.0957 

6 

63 

.2743 

14 

4 

79 

1.0920 

5 

71 

1.1109 

6 

65 

.3230 

ai4 

4 

79 

1.0920 

5 

71 

1.1109 

5 

73 

.8793 

16 

3 

67 

1.0573 

5 

77 

1.1208 

5 

59 

.6914 

18 

3 

71 

1.0558 

5 

79 

1.1215 

5 

63 

.7597 

^18 

3 

71 

1.0558 

5 

79 

1.1215 

4 

67 

.9717 

20 

3 

75 

1.0536 

4 

71 

1.0758 

5 

65 

.7889 

22 

3 

79 

1.0508 

4 

73 

1.0681 

5 

67 

.8153 

^22 

3 

79 

1.0508 

4 

73 

1.0681 

4 

75 

.9862 

24 

3 

101 

1.0346 

4 

75 

1.0606 

5 

71 

.8602 


^Relation (40) was used in addition to relation (37). 
^Unchanged when using relation (34) rather than relation (35). 


poorer "necessary but not sufficient" condition. The net effect, as seen in table I, is 
generally to keep N larger and m smaller than N and m for the other wings, for 
a fixed value of This amounts to putting the first chordwise control point nearer the 

leading edge while slightly increasing the spanwise distance between tipmost integrated 
and integrating stations. 


Another procedure was implemented which extended relation (36) to the following 
more stringent form: 

lY. 


^0 > 


tan A 


le 


(40) 
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This relation biases Xq away from the leading edge by weightii^ the right-hand side of 
the inequality. Relation (40) is used in addition to the "necessary but not sufficient" rela- 
tion (37). Relations (40) and (37) were used for the following values of Cj^: 10, 14, 18, 

and 22; these resvilts are identified in table I. Using relations (40) and (37) had no effect 
on the wing 1 and 2 solutions; however, there were improvements for wing 3. Note that 
the Cj-j i values more quickly approach unity with a smaller product of N 

and m by using relationships (40) and (37) than with relationships (36) and (37). 

Effect of Initial Values of N and m 

Thus far the results have shown that (1) Cj « 10 is a good compromise choice 
for untapered wings regardless of the relationships used to determine satisfactory 
N,m sets and (2) for tapered wings, relations (40) and (37) yield improved results. Con- 
sequently, a hypothesis is offered that the N,m sets which result from all preceding 
relationships may be dependent on the initial values of N and m. This dependence 
is suspected because even with local satisfaction of the relationships, there were 
Cd ii/Cj) values in table I far from unity for raises of Cj. It appears that both 
"local" and "global" conditions must be satisfied. Therefore, depending on the initial 
values of N and m, or where the procedure starts from, both conditions may not be 
satisfied. 

To examine this hypothesis, a study was conducted based on relations (40) and (37) 
with Cj = 10. Results of this study are summarized in table n for 20 sets of initial N 
and m values. 

Before the results are analyzed, the computerized procedure for determining the 
N,m sets should be explained. It is as follows: 

(1) The initial guesses for N and m are used in relation (40), in conjunction 
with equations (38) and (39); if the relation is not satisfied, m is increased by 4 and 
another attempt is made. 

(2) Increasing m by 4 continues until relation (40) is satisfied or m > 101 in 
which case m is set to 101, N is reduced by 1, and the code attempts to satisfy rela- 
tion (40) again. 

(3) If this reduction in N is not successful in satisfying relation (40), then N is 
again reduced by 1, still with m = 101, and the process is repeated until the relation is 
satisfied or N = 0. If N = 0, the program writes a message and stops. 

(4) Once relation (40) is satisfied, the N,m combination that led to satisfaction of 
relation (40) is used in relation (37). If either the N or m is not appropriate for 
relation (37), it is incremented as described previously, and the procedure returns to 
step (1) where relation (40) must once again be satisfied. 
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TABLE n.- EFFECT OF INITIAL N AND m VALUES 
jCj = 10; M = 0 


Initial value 


Wing 1 


Wing 2 


Wing 3 

N 

m 

N 

m 

^D,ii/^D,i 

N 

m 

^D,ii/^D,i 

N 

m 


2 

11 

2 

39 

1.1108 

2 

31 

1.0271 

2 

31 

0.9341 

2 

31 

2 

39 

1.1108 

2 

31 

1.0271 

2 

31 

.9341 

2 

51 

2 

51 

1.0863 

2 

51 

1.0057 

2 

51 

.9401 

2 

71 

2 

71 

1.0607 

2 

71 

1.0119 

2 

71 

1.0047 

2 

91 

2 

91 

1.0501 

2 

91 

1.0177 

2 

91 

1.1185 

4 

11 

4 

67 

1.0746 

4 

51 

1.1188 

4 

51 

.8840 

4 

31 

4 

67 

1.0746 

4 

51 

1.1188 

4 

51 

.8840 

4 

51 

4 

67 

1.0746 

4 

51 

1.1188 

4 

51 

.8840 

4 

71 

4 

71 

1.0827 

4 

71 

1.0758 

4 

71 

.9804 

4 

91 

4 

91 

1.0939 

4 

91 

1.0097 

4 

91 

.9858 

6 

11 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

6 

31 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

6 

51 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

6 

71 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

6 

91 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

8 

11 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

8 

31 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

8 

51 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

8 

71 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 

8 

91 

4 

67 

1.0746 

5 

59 

1.0380 

5 

63 

.7597 


(5) After both relations are satisfied, a test is made to determine whether the prod- 
uct of N and (m + l)/2 is larger than 200, the maximum allowed. This product is the 
unique number of modal factors on the wing under either symmetrical or antisymmetrical 
loading conditions. 


(6) If the product is less than or equal to 200, the aerodynamic solution is started. 
If it is greater than 200, N is reduced by 1, m is reduced to 


2 


fc 


I 


Greatest whole integer 



- 1 


and the procedure returns to step (1) for another start. 
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Examination of table II prompts the following observations: 

(1) A few initial N and m values satisfied both relations and became the final 
N,m sets. 

(2) For N = 6 or 8 and any m, the same final N and m values for a particular 
wing were obtained. 

(3) The best overall Cp ^ results were obtained with initial values of N = 2 

and m = 71, and second best witn initial values of N = 4 and m = 91. 

(4) The delta wing (wing 3) shows the most sensitivity in Cjj . to the initial 

N and m values. 

(5) For N = 2, Cjj i convergence generally occurs as initial values of m 

increase up to 71 for all wings and up to 91 for wing 1. For N = 4, Cjj njCjy . con- 
vergence occurs as initial values of m increase up to 91 for wings 2 and 3, whereas a 
mild divergence occurs for wing 1. 

Computer Requirements 

An overlayed computer program which permits up to 200 flat -wing modal factors 
and m = 101 requires approximately 77000g words in central memory for the sym- 
metric mode of operation and more in the antisymmetric mode. These computer require- 
ments for a potential-flow aerodynamic solution were considered excessive; thus, an 
operational version was developed which permits up to 100 flat -wing modal factors, 
m = 41, and requires only 510000 words of central memory in either mode of operation. 
The consequences of using this computer program are examined before the results 
obtained with it are discussed in the remainder of this report. 

It should be pointed out that the procedure for finding a final N,m set from the 
initial values is slightly different in the smaller program. The difference is primarily 
that a value of m considered "large enough" (usually 41) is specified initially, so that 
only reduction in N is permitted. In general, this is satisfactory as is seen 
subsequently. 

The primary consequence of using the smaller program is that for some wings a 
value of N which is too small results from relations (40) and (37). This leads to the 
boundary condition being satisfied only at as few as one chordwise position. The only 
way that this situation could be remedied is to increase the value of m allowable, 
essentially to return to the original program discussed. The operational computer pro- 
gram, with modest computational requirements, is capable of generating acceptable 
Cq ii/^D i results for most wings, whereas the original computer program with larger 
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computational requirements would provide acceptable Cj^ . results for an unknown 

number of additional wings. Thus, the smaller program was chosen. 

Comparison With Previous N,m Set Criterion 

Reference 7 presents an approximate formula for the relationship between m 
and N which can be used to obtain results that are considered converged. This relation- 
ship was developed for rectangular wings and is 

m » (4 to 5) |3A ^ (41) 

where A is aspect ratio. It is of interest to compare the formulation of reference 7 
with relation (33) for the same wings. For rectangular wings, relation (33) reduces to 

Xu » |Yo| (42) 

By using equations (38) and (39), relation (42) can be rewritten as 


1 - 



> C. 


/3 /b' 

2 c(y)\2 



which with reduction becomes 


1 - cos 


2t7 


2N + 1 


>C. ^ 


cos ( — - cos \ 

\m + 1/ \ m + 1/^ 


(43) 


(44) 


where the aspect ratio A is 

A = ^ 

2 c(y) 


For Qn/{m + 1^ < 1, the right-hand side of relation (44) can be expanded to yield 


1 - cos 




sr 


2(m + 1)2 


(45) 
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and so, 


m> — 


Cj l3A37r‘‘ 


1 - cos 


2n 


2N + 1 


- 1 


(46) 


As an example of the different results obtained from these two procedures, consider the 
following values for the variables: ;S = 1, A = 2, N = 4, and = 10. Then equa- 

tion (41) yields 


m « 10 


and equation (46) yields 
m > 34 


ADDITIONAL NUMERICAL STUDIES 

Two additional numerical studies were made. The first examines the effect of the 
new LSC on the ratio of near -field to far-field vortex drag, and the second examines the 
effect of on the overall longitudinal aerodynamic results. 

Effect of LSC on Vortex Drag Ratio 

To show the improvement obtained by using the new LSC, the vortex drag ratios for 
two planforms, rectangular and delta, are studied over a range of N and m. Results 
using the LSC described in reference 7 are compared with those using the new LSC. It 
should be stated at the outset that for either LSC, vortex drag ratios near 1 are not to be 
expected for any arbitrary combination of N and m. Figure 1 presents the results of 
this numerical experimentation. The general conclusions for both planforms are that the 
new LSC (1) provides a stable solution which converges for a fixed value of N with 
increasing m and (2) yields convergence with increasing N only for a sufficiently 
large value of m. 

The second conclusion is strongly related to relation (37) and signifies how viola- 
tions of the N,m criterion decrease the solution accuracy. For example, if N is too 
small in relation to m, the downwash equation cannot be satisfied more accurately in the 
chordwise direction. If N is large, the misrepresentation of the logarithmic singularity 
near the leading edge reduces the accuracy, depending on the seriousness of the violation. 
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On the other hand, larger N can better satisfy the boundary condition so that the accu- 
racy is increased. The whole problem is the following; N and m shovild be reason- 
ably large so that the boundary condition can be better satisfied, while simultaneously the 
correct representation of the logarithmic singularity should be preserved. The N,m 
criterion (relation (37)) deals with the latter condition only. Table II illustrates the point 
that solutions for small N may be as good or better than those for larger N where 
each has an m value satisfying relations (37) and (40). 

Effect of on Longitudinal Aerodynamic Results 

Results of the second numerical study are presented in table m. In it the aerody- 
namic results obtained by employing both the old LSC from reference 7 and the new LSC 
for a variety of planforms are given. Reference 7 used the N,m relationship given by 
equation (41). Results using the new LSC are presented for two values, and rela- 
tions (40) and (37) were employed for the N,m relationship. 

Upon comparing the results in table HI, the following observations can be made: 

(1) The new values of are generally greater than those of reference 7. 

(2) The Cj^ values vary only slightly with solution or Cj . 

(3) The = 10 solutions generally have N values greater than the = 14 
solutions. 

(4) The Cp I values for Cj = 10 are often slightly farther from 1 than 

those for Cj = 14, with the exception of the more complex planforms. 

(5) At N,m combinations based on equation (41) the drag ratios for the new LSC 
solutions are generally farther from 1 than those obtained with the N,m combinations 
based on relations (40) and (37). 

These observations indicate that the Cj = 10 solutions are almost as good as those for 
C]^ = 14 and generally provide one more chordwise point for local downwash satisfaction. 
Hence, Cj = 10 is used subsequently. 

It should be noted again at this point that with m increased the new version would 
produce improved results. This is not necessarily true with the old version, as has been 
shown in figure 1. 



TABLE m.- LONGITUDINAL AERODYNAMIC RESULTS USING OLD AND NEW LSC FOR SEVERAL PLANFORMS 


Wing type 

Aspect 

ratio 

A, 

deg 

X 

M 

Cl 

N 

m 

C 



Cd.U 


Cd.u/Cd,! 

new LSC and 
equation (41) 

Old LSC 
(a) 

New LSC 

Old LSC 
(a) 

New LSC 

Old LSC 
(a) 

New LSC 

Rectangular 

2 

0 

1.0 

0 


4 

11 

0.0432 


0.208 


0.9943 


0.9549 






10 

4 

41 


0.0432 


0.210 


1.0337 







14 

3 

41 


.0431 


.210 


1.0288 


Rectangular 

7 

0 

1.0 

0 


4 

37 

0.0768 


0.239 


1.0277 


0.8414 






10 

2 

41 


0.0777 


0.244 


1.1067 







14 

1 

41 


.0764 


.250 


1.3209 


Sweptback 

2 

45 

1.0 

0 


4 

11 

0.0398 


0.170 


0.5311 


0.5371 






10 

4 

41 


0.0398 


0.171 


1.0610 







14 

3 

41 


.0396 


.169 


1.0297 


Sweptback 

6 

46.17 

0.60 

0 


4 

31 

0.0615 


0.252 


-0.5252 


-0.5573 






10 

2 

41 


0.0607 


0.268 


1.0818 







14 

1 

41 


.0596 


.263 


1.0363 


Delta 

2 

63.4 

0 

0.13 


4 

11 

0.0391 


0.363 


-0.8669 


-1.5331 






10 

3 

41 


0.0385 


0.390 


0.9721 







14 

2 

41 


.0385 


.389 


.9372 


Cropped 

2 

45 

0.25 

0 


4 

9 

0.0432 


0.248 


0.5011 


0.3748 

diamond 





10 

4 

41 


0.0430 


0.266 


1.0749 







14 

3 

41 


.0430 


.265 


1.0009 


Backward, 

2 

0 

0.25 

0 


4 

11 

0.0426 


0.191 


0.9567 


0.9813 

cropped delta 





10 

4 

41 


0.0424 


0.193 


1.0325 







14 

3 

41 


.0423 


.191 


1.0355 


Double delta 

1.415 

82/60 

0.10 

0 


4 

9 

0.0317 


0.433 


-0.4254 


-0.8222 






10 

4 

41 


0.0313 


0.459 


1.1348 







14 

3 

41 


.0313 


.461 


1.0706 


Double delta 

1.7 

75/60 

0.10 

0 


6 

13 

0.0365 


0.375 


-0.9101 


-1.8484 






10 

3 

41 


0.0362 


0.388 


1.0462 







14 

3 

41 


.0362 


.388 


1.0462 


Double delta 

1.97 

. 45/60 

0.10 

0 


6 

15 

0.0407 


0.324 

1 

-0.2372 


-1.2064 






. 10 

3 

|41 , 


0.0404 


0.334 


1.0735 






1 

1 

2 

|41 1 


.0405 

i 

.335 , 


.9196 



oi ^aken from reference 7. 
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APPLICATIONS 

Table IV presents the vortex drag ratios for a wide variety of planforms and sub- 
sonic Mach numbers. All results are obtained from the hands-off computer solution, 


TABLK IV.- VORTEX DRAG RATIOS FOR A VARIETY OF PLANFORMS 

AND MACH NUMBER 

[Ci = l^ 


Wing type 

Aspect 

ratio 

A, 

deg 

X 

M 

N 

m 


Rectangular 

0.2 

0 

1.0 

0 

8 

23 

0.9799 

Rectangular 

.3 

0 

1.0 

0 

7 

23 

1.0956 

Rectangular 

.4 

0 

1.0 

0 

6 

23 

.9721 

Rectangular 

1.0 

0 

1.0 

0 

5 

31 

1.0344 

Rectangular 

2.0 

0 

1.0 

0 

3 

31 

1.0350 

Rectangular 

3.0 

0 

1.0 

0 

3 

41 

1.0329 

Rectangular 

7.0 

0 

1.0 

0 

2 

41 

1.1067 

Sheared rectangular 

3.5 

0 

1.0 

.3 

3 

41 

1.0366 

Sheared rectangular 

3.5 

20 

1.0 

.3 

3 

41 

1.0503 

Sheared rectangular 

3.5 

40 

1.0 

.3 

3 

41 

1.0734 

Sheared rectangular 

3.5 

50 

1.0 

.3 

3 

41 

1.1348 

Sheared rectangular 

3.5 

60 

1.0 

.3 

2 

41 

.9751 

Sheared rectangular 

3.5 

70 

1.0 

.3 

1 

41 

1.1589 

Sheared rectangular 

3.5 

75 

1.0 

.3 

1 

41 

1.1241 

Sweptback 

1.0 

45 

1.0 

0 

4 

41 

1.0179 

Sweptback 

2.0 

j 

45 

1.0 

0 

4 

41 

1.0610 

Delta 

4.0 

45 

0 


3 

41 

.9705 

Cropped delta 

3.273 

: 45 

.1 

.6 

3 

41 

1.0764 

Cropped delta 

2.667 

: 45 

•2 

.6 

4 

41 

1.1397 

Cropped delta 

2.154 

45 

.3 

.6 

4 

41 

1.0648 

Cropped delta 

1.714 

45 

.4 

.6 

4 

41 

1.0279 

Cropped delta 

1.333 

45 

.5 

.6 

4 

41 

1.0197 

Cropped delta 

1.668 

63 

.1 

.6 

3 

41 

1.0378 

Cropped delta 

1.359 

63 

.2 

.6 

4 

41 

1.0986 

Cropped delta 

1.097 

63 

.3 

.6 

4 

41 

1.0207 

Cropped delta 

.873 

63 

.4 

.6 

4 

41 

.9887 

Cropped delta 

.873 

63 

.4 

0 

4 

41 

1.0313 

Cropped arrow 

1.069 

63 

.54 

0 

4 

41 

1.0148 

Cropped arrow 

1.917 

63 

.29 

0 

3 

41 

1.0165 

Cropped diamond 

.738 

63 

.32 

0 

4 

41 

1.0248 

Arrow 

3.25 

60 

0 

0 

2 

41 

.9397 

Diamond 

1.75 

60 

0 

0 

3 

41 

1.0379 


26 



Ill 


with the N,m sets derived from relations (40) and (37) and equations (38) and (39) with 
Cj = 10. Many of the initial N and m values used in the solutions were N = 4 and 
m = 41. From table IV, it can be seen that generally the near -field drag is within 
±10 percent of the far -field drag, and in many cases the error is even smaller, within 
±5 percent. The larger errors do not occur in any systematic pattern that could be asso- 
ciated with a particular wing type, aspect ratio, sweep, taper ratio, or Mach number, nor 
do the values of N and m selected by the program lead to any correlative trends. It 
is expected that these errors can be reduced if m was increased beyond the 41 limit 
imposed by the smaller computer code. 

As noted in the discussion of relation (33), the computer program may have diffi- 
culty in satisfying that relation for wings with higher leading-edge sweep. An example of 
this occurs for the sheared wing results given in table IV. A solution for A = 80° was 
attempted but the conditions in relation (33) could not be satisfied for m = 41 or 
m < 41 at any N value. Hence, it was not included. This led to the use of A = 75° 
which was not too large to satisfy the aspect ratio, taper ratio, and Mach number require- 
ments inherent in relation (33). 


CONCLUDING REMARKS 

This report has presented the development of a new logarithmic -singularity correc- 
tion term to be used in Multhopp lifting-surface computer programs. One novel aspect 
of the correction term is that it is expressed as a function of the complementary modulus 
of the complete elliptic integrals. In this form the correction term contains both chord- 
wise and spanwise variations in the integrand near the spanwise singularity. It also 
leads to the establishment of a relationship between the chordwise control point nearest 
the wing leading edge and the spanwise distance between integrating stations and control 
points, so that the correction term computed will be valid. With this term set, reliable 
aerodynamic results can be obtained in a single computer pass. In addition, even if the 
previously mentioned relationship is not used, it has been determined that the new correc- 
tion term leads to stable converging solutions for a fixed number of chordwise control 
points with increasing number of spanwise integration stations. Stable converged solu- 
tions have not always occurred for arbitrary planforms and subsonic Mach numbers. 

Numerical studies indicate that with this new relationship, the ratios of near-field 
to far -field vortex drag for a variety of planforms and subsonic Mach numbers are gen- 
erally between 0.90 and 1,10 and in many instances between 0,95 and 1,05. 

Lai^ley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 15, 1977 
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APPENDIX A 


EVALUATION OF J^, Jg, AND dk’/d77 
Evaluation of 

The closed-form expression for has been obtained by Wegener and published 
in reference 10. However, details of the derivation are not available. They are included 
here for the sake of completeness. In this appendix, all page numbers refer to 
reference 11. 

can be written as 



V ^1 ‘ “ ^Xi(l -Xi)(/(X-Xi)2+Y= 


dXi 


-(1+x) 

-^0 


pi Xj dXj p 

.1 Xj^ dXj 

“ \/Xi(1-Xi)^(X-X,)2+y2 ' 

^ \/5^('-Xi)^(X-X,)2*y2 

+ Jj3 + J^4 

(Al) 


From p. 133, some parameters involved in the results are defined as follows: 

A^ = (1 - xf + (A2a) 

(A2b) 
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APPENDIX A 


, 2 _ 1 - (A - B)^ 
4AB 


(A2c) 


g = 


\Jab 


(A2d) 


(p - TT 


(A2e) 


The first integral in equation (Al) can be shown to be tt/2. Now, by Item 259.00 

on p. 133, 


'12 


= x 

*^0 


dX, 


^Xi(1-Xj)^(X-Xj)2+y2 


= Xg F(<p,k) = 2Xg K(k) 


(A3) 


Again, from Item 259.03 on p. 133, it is found that 


■’13 ' 


Xj dXj^ 


)1Xi(1-Xi)^(X-Xi)2..y2 


= -g(l + X) 


n 


2B n( ,k) 


A - B A - B 


(A4) 


where f j, defined in Item 361.54 on p. 215, equals zero and 

A 

«2= 


a = 


A - B 
A + B 


l-«2= 4AB 


(A + B)^ 


0-2 _ (A - B)^ 


«2 _ 1 


4AB 


(A5) 
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Similarly, 


'14 


*^0 


Xi^ dXj 




gB^ 


F((p,k) + 2fl!2(a - ffg) ^ ^ ' " 2 )^ ^2 


(A - B)2 

where from Item 341.04 on p. 206, 

Q!^(2k^ - l) - 2k^ 


R2 - 


(a^ - l)(k2 + a\'^) 


R, + 2k2R . - k^R „ + ^ 

■*■ 1 + o! cn u 


I 2 . 2 

where cn u = cos cp, sn u = sin cp, dn n = \J1 - k sin <p, and 


(p=n 


Ri =-^^uicp,-^ — ,k) 


" 1 - 0(2 
R_1 = F(<p,k) 


«2_i 


R 9 = ^ 
-2 k2 


(k2 - Q!2k'2) F(^,k) + q' 2 E((/J,k) 


Hence, 


= 


__gB‘ 


(A - B)2 


I 1 - a‘^ \ a“^ - 1 




{“ 

-“ 2 )^ 

«2| 

( 2 k 2 - 1 ) - 2 k^ j 3 | 


{J - 1 ) 

(k2 + Q!2k’2) 

._ 

l-«2 

iS^? 0 

^ «2 _ 1 y 


2k2 F((/j,k) - (k2 - a^k'^) F((^,k) - E((p,k) 
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By substituting equations (A2) and (A5), becomes 


Jl4 - g 


K(k) + 2AB E(k) + - ^ ^ [2 - (a^ - B^ 

A - B 2 A - bC ' 



(A6) 


Combining the expressions for Jj^ 2 > Ji 3 > ^14 (A3), (A4), and (A6)) 

gives 


J + 2Xg K(k) - g(l + X) 

Lt 


2B 

A - B 


K(k) + 


A + B 
A - B 



K(k) + 2AB E(k) + u( 

A - B ' A - B \q,2 _ 1 ’ y 


-|(A.B)2 



= I + 2g K(k) + AB E(k) 

2 A - B 


1 A + B 

2 A - B 



- j(A + B)2 

4 



From Item 410.01 on p. 225, 



can be written as 


n 



= 4ABk^ K(k) + 


o l (A - B)^ (A + B)^ 1 


4AB 


4AB 4AB 


= 4ABk^ K(k) + ^\/^|a - b [ a (i//»,k) 
|A + B] 0'*' ’ ' 


(A7) 


where 


\p' = sin 


^|a - b] 
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By equation (A7) , Jj becomes 


j = I + K(k) + AB E(k) - I X 

1 2 ^A - B 2 A - B 


4ABk^ K(k) + tt>/aB 


A - B 


A + B 




-^(A + B)2 
4 


4ABk^ K(k) + tt^AB 



- 


A-B 

AQir,k) 

\ 

A + B 

J 


|.2g 


K(k) + AB E(k) - AB(A , + B) _ WabJA_^ Ajxly' ,k) 

A - B A - B 4(A - B) ” 


(A8) 


It is known (see p. 36) that for any ijy' , 

AQ(-)//',k) = -Ao(t//’,k) 

Hence, 


A - B 
A - B 


A^ip^k) = AQ(i//,k) 


where 


\jy = sin"^ (A - B) 


The coefficient of K(k) in equation (A8) can be simplified as follows: 

XA _ AB(A + B) jj2 ^ XA _ AB(A + B) 1 - (A - B)^ 
A-B~ A-B A-B A-B 4AB 


(A9) 


2a[i - (a^ - B^)) - (A + B) + (a2 - B^) (A - B) 
4(A - B) 

^ -a(a^ - B^) - b (a^ _+ A - B 

4(A - B) 


= z(A^tB)i±l_ABk’2 
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Thus, substituting for g from equation (A2) reduces to 
Jj = I + 2)/AB(E(k) - k'2 K(k3 - I A(j(;^,k) 


= I + i\/AB(E(k) - k’2 KCkjj 


(AlO) 


Evaluation of J 2 


Jg can be written as 


•^0 ^X,(l-Xi 


X - Xi 


=1 


1 dX 


^(X - Xi) 

i= + X 


2 + 


dX, 


dXi 






X^ dX^ 




The first integral in the expression for J 2 can be shown to be n. The last two inte- 
grals can be expressed in the same way that J ^2 a-iid Jj^g were (eqs. (A3) and (A4)). 
The resulting expression for J 2 is 


J 2 = 7T + 2Xg K(k) - g 


K(k) + 

A - B 


A + B 
A - B 



Using equation (A7) for 



results in 


J 2 = 7T + 2Xg K(k) - g 


K(k) + ^ + ® 


A - B 


B 


4ABk^ K(k) + TT yjA^ | ^ ~ 

lA + Bl 
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The coefficient of K(k) can be reduced to 


2g 


^ B 2ABk^(A + B) 
ATi' a-B 


= 2d 


2XA - 2XB + 2B - (A + B) + (a^ - B^)(A - B) 
2(A - B) 


= 2g 


2XA - 2XB - 2B - A - B + (1 - 2X)(A - B) 
2(A - B) 


= 0 


It follows that 

Jg = wQ - 7VQ(i//,k)] (All) 

where xj/ is defined in equation (A9). 

Evaluation of dk'/dq 

For 77 — y ^i.e., X » Xq^, the following relation can be obtained from the defini- 
tions of X and 0: 


i(l 

2 


cos (ft) = 


^ - ^le^y) 

2c(y) 


Hence, x is related to (f) through 
X = lig(y) + c(y)(l - cos (^) 

Thus, from the definitions of X and Y, 

_ 1 ^le^y^ ■ ^le^’^^ 

^ 2 c(r/) 

Y - 1 j 3 (y - 77) 

2 0(77) 


(A12) 


(A13) 
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Differentiation of equation (A13) with respect to rj gives 


[ileW ^ 

cHv) 


iv) 


dri 2 c(? 7 ) 


drj 2 


1 _ (y - v) c'jri) 
civ) c 2 ( t 7 ) 


(A14) 


where c'(? 7 ) and ^iQiv) can be expressed in terms of the sweep angles as follows: 
c'iv) = |(tan A^g(? 7 ) - tan Aig(ri)) 




(A15) 


From definition, 


k’2 = 1 


,2 _ (A + B)^ - 1 
4AB 


- X + ^ 1 

2AB 2 


(A16) 


Now, differentiation of equation (A16) gives 


2k’ 


dk» 

dr] 


ov dX ov dY / o o' 

d^ ' d;r d^ . (x - X + y2. 


2AB 


2A^b2 


B 


But from the definitions of A and B 


dA^^^\ 

dr] drj ) 


2A — = -2(1 - X) — + 2Y — 
d ?7 drj dr] 


2B ^ = 2X -— + 2Y — 


dt] 


dr] 


dn 
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dj] dr] AB \ drj 


— Cb^(1 - X) — + B^Y — + A^X — + A^Y — 
^BL <^’7 *7 


^BL -J drj AB dt] 


Also, 


X^ - X + Y^ = AB(2k’2 - l) 


By these relations, it follows that 


, dk’ _ 2X - 1 dX . Y dY 


dr? 2AB dq AB drj 


AB(2k'^ - l) / 1 
2A^b2 1 AB 


:a2+b2 


')x - B^ 


^ + Y — 
dt] AB dt] 


The coefficient of dX/dq in equation (A17) can be written as 


2X - 
2AB 


_ j (2k'2 - i)[(a^ + B^)x - B^ 


2A^b2 


= _^^(2AB + A^ + B^) - AB - B^ - 2k’2(A^ + b^)x + 2k'^B^ 
= _^^(4ABk’2 + l) - AB - b 2 - 2k'2(A^ + B^)x + 2k’2B^ 

= — L(-2Xk'2(A - B)^ + 2B^2 - 2ABk’2 

2a2b^^ L -* 
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Similarly, the coefficient of Y — in equation (A17) can be written as 

dT? 


1 (2k»^ - i) (a^ + B^) _ 1 

AB 2A^b2 2a2r^ 


2A 


j^A + B)^ - 2k’2(A^ + B^^ 

- 2k'2(A - B)2] 


Substitution of these results into equation (A17) gives 


2k* ^ = 


k’^ 


^7 A^B^ 


X(A - B)^ + B^ 



d?7 


+ 


2A^b2 




- 2k’^(A - B) 


3 


;3(y - 7 ) dY 
2 c(? 7 ) drj 


(A18) 
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DERIVATION OF LSC 


Before making any expansion of Jj and Jg for small Y, it is important to know 
the relation between and k’^ as defined in equation (A16). From equation (A16), 

1^,2 . - X + .1 

2\J(1 - xf + Y^ /x^ + Y^ ^ 

2 

Note that k'^ approaches zero when Y vanishes. This expression can be reduced to 
a quadratic equation in Y^ as follows: 





+ (1 - X)^ 



Y 


2 


4k’2(k’^ - l)x2(l - X)2 = 0 


(Bl) 


This algebraic equation can be solved for exactly and then the results expanded for 

small k’^ or it can be solved by the perturbation method by noting that both Y^ 

2 

and k’ are small. By using the latter method, it is observed that to the first 
approximation, 

Y^ = 4k'^X^(l - X)^ 

Since the first term in equation (Bl) is o(k'®), it may be neglected. Thus, by retaining 
terms involving k'^. 


Y^ = 


4k’^(k’^ - l) X^(l -X)^ , o(l:’6) 

1 + 4k’2(k'2 - l)jx2 + (1 - X)^ 


= 4k'2x^(l - X)^ 1^1 + k'^d - 8X(1 - X^^ + o(k'®) 


(B2) 


From equation (B2), it is seen that small k' implies always small Y, but the converse 
is not true. Hence, it is more natural to develop Jj^ and Jg for small k' rather 
than for small Y. 


38 



APPENDIX B 


It is known that (see p. 36 of ref. 11) 


Ao(i//,k) = |(E{k) F(i^,k’) + K(k) E(i//,k’) - K(k) F(i//,k')j 


(B3) 


From pp. 299-300 of reference 11, the following expansions for small k' are valid: 

E(k) = 1 + In -1 + . . . 

2 k' 


K(k) = lnA^l+^ + . . .j H-. . . 

F(i//,k’) = i// + i k’^ (i// - sin i// cos il/) + . . . 
4 

E(;//,k') = i// - i k'^ (i// - sin \p cos \}/) + . . . 


(B4) 


where only the logarithmic components of E(k) and K(k) have been retained. Substi- 
tution of equation (B4) into equation (B3) gives 


U ^ 



^ " si" 



k*' 


\p — “ sin i// cos \p) - Ip 


k’ 


^ - sin \p cos 


Therefore, 


Logarithmic 

component 


of Aq 


^^-sin 4/ cos ip In k’^ + o(k'^ In k') 


From equation (A9), sin = A - B, so that 


Logarithmic •, i 9 « / 4 \ 

component = -—(A - B)\Jl - (A - B)^ k’^ In k' + o(k*’ In k'j (B5) 

of Aq ^ 
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It follows that (see eq. (16)) 


Logarithmic / s „ / a \ 

component = (A - B)\/l - (A - B)^ k*^ In k' + o(k*'* In k'j 
of Jg 


Similarly, by using equation (B4), equation (15) becomes 


J 

'^ 1~2 


1 - A(,(v-,k) f I mil * V F ■ 


(B6) 


Then, 


Logarithmic 
component 
of Jj 


iji (A - B)\Jl - (A - B)^ k'2 In k’ - | \/ab k’^ ^ + o(k'‘^ In k’) 


ij^A - (A - B)^ + 2\[aB k’^ In k’ + o(k’^ In k’) 


(B7) 


Equations (B6) and (B7) can be further expanded by noting that 


A = J(1 - X)^ + = 1^1 - X)^ + 4k’^X^(l - X)^ + o(k’^) 

= 1 - X + 2k’^X^(l - X) + o(k*'^) 


1/2 



= X + 2k’^X(l - X)^ + o(k’^) 


/ (B8) 


where equation (B2) has been used for y 2. Thus, by retaining only terms of the order 
of k*^ In k’, equations (B6) and (B7) become 
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Logarithmic , , ^2 / “i \ 

component = 2 JXq(1 - Xq) (1 - 2Xq) In k’ + o(k’'^ In k'j 


of J, 


= - ^o) (1 - 2Xo) 1” ■‘■((y - ■»' 


Logarithmic , , ^ q \ 

component = 2 JXq^1-X^(1 - Xq^ k In k* + o(k’ In k'J 


of J 




where Xq, the value of X at v = Y, has been used to replace X. The derivative 
dk'/dT? is evaluated in appendix A. It is known from equations (A14) and (A18) that 

^ &-IJ21 + o(k') 

4k’A'^B^ 4 cHv) 


From equation (B2), 


y2 = ^(y - ^) - = 4k’^X^(l - X)^ 1 + 
4 c2(tj) *- 


o(k'^)' 


y _ 77 = s ign (y - 7?) 4 c^(77) _ k^X(l - X ) ^ o(k’2)J 

4 c(y) k*XQ(l - Xq)[i + o(k’2^ 


Hence, 
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dq 

dk’ 


- ^o) 


_4 c^^(y) ^ 


)3^4 c(y) k’XQ^l - Xq) sign (y - t?) jl + o(k’^^ 


+ 0(k') 


= . signJyjL^Z) 4 c(y) Xq(1 - Xq) + 0(k’) 


(B9) 


Furthermore, 


^ !— ,= -’!) 4c(y)X()(l -X(|) 

<a=’ (y - <if ^ ' 


16 c2(y) k'2x„2(l - X„)2 



= - sign (y - q) 


^ 

4 c(y) k’2xQ(l - Xq) 



It follows that 


Logarithmic , ,9 9 Htt 

component = 2^q(1 - Xg) (1 - 2Xq) k’'* In k’ (y - riy ^ 


-sign (y - n) 


4 c(y) k’2Xo(l - Xq) 


*°(f) 


■kO(k'^ Ink') 


= -sign (y - 7j) |3 (y - vf f^lnk’ + o(k’2 In k’) 

2 c(y)^Xo(l-Xo) ’’ 

= -sign (y - 7 ?) (y . i ,)2 ^ In k- -k o(k'3 in k’) 

c(y) sin 0 d?? 


(BIO) 


Similarly, 
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Logarithmic 


= 2(1 - Xq) - Xq) k-2 In k' ^ (y - 77) 


-sign (y - 7j) 


4 c(y) k’2xQ (1 - Xq) 


+ o(k’^ Ink’) 


= -sign 




-sign (y - rj) — 2_ ^ <p ^ j^, ^ q(i ^.3 j^,) 

2 c(y) sin 0 dr; 


If equations (BIO) and (Bll) are substituted into equation (14), the logarithmic 
ponent of Ij can be shown to be 




g^(0) + (1 - X) g2(0) 


-ign (y . ,) -2^ (y . ,)2 ^ in k’ H. o(k’3 In k-^ 


^^^Csign (y - T?) ^ ^ (y - v)^ ^ + o(k'^ In k')| 

L 2 c(y) sm cp dq J 



(y . 77)2 f: In k- 

= -sign (y - ,) I -y-2 Gj(7,) 


where 


G.(<^) = j0 + cos (j + l)^ ^ ^ j ](p + a +1) sin (j + 1)0 

J sin 0 

It should be noted that equation (B12) can be reduced to the conventional form, 
as that used by Wagner (ref. 5), by' simply using equation (B9) for dk’/dr]: 


(Bll) 

com- 


(B12) 


(B13) 

such 
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^ = -sign (y - v) 
dr] 


|8 


4 c(y) Xq(1 - Xq) 


= -sign (y - v) — — ^ « 
c(y) sin^0 


and to the first approximation, 


In k’ = In 


4x2(1 - X)^ 

= ln |Yl-ln[2Xo(l -Xq) 


If only the term associated with In | Y| is retained, it is seen from equation (Bl 2 ) that 


Lj = -sign 


(y - J?) I In I Y| f-sign (y - 77) ^ - "1 GM) 

^ c(y) [_ c(y) sin2<;^ •> 


= i 1 Y |2 In |Y| + o(|Y |3 In |y|) 

^ sin^(p ' 

where c(y) has been expressed in terms of c(r]) through the relation 
c(rj) = c(y) + c’(y) (?7 - y) + 0(?? - y)^ 


Since 


dh 


1(^1) 


dXi 


01=0 


^ of cos + cos (j + 1)0 ^ ^ j sin j0 + (j + 1) sin (3 + 1)0 | 

L sin^^ sin^(p J 


o 

■ 2-2 

sin'^0 


it follows that 




dXi 


1 y |2 In |Y] 


01=0 
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which is the same expression as that used by Wagner (eq. (16b) of ref. 5) except for the 
factor 2 /jt which is due to the different definition of hj used here in this report. 

To perform the spanwise integration defined in equation (5), the chordwise influence 
function Ij(X,Y,t?) and the loading function qj(Tj) Ij(X,Y,? 7 ) may be decomposed as 
follows: 


I.(X,Y,7?) = l^*{X,Y,v) + L. (B14) 

q.iv) Ij(X,Y, 77) = q.(y) (l. - I.*^ + j~q.(77) - q.(y) I. + q.(y) I.* (Bl5) 


When equations (B14) and (B15) are substituted into the spanwise integral in equation (5), 
the first term in equation (Bl5) can be exactly integrated and gives 



The notation has been made compatible with Multhopp's notation (ref. 1). Reference 7 
gives a detailed explanation of the subscripts which have been introduced. It follows that 


-2vrw(0g,yJ 


part 1 



3=0 


277 c^yj 


(‘^s) (^i In 



(B16) 
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where 




,.,2 

1 


.j2=(l-Xlf .V 


* y/ 


^l = t 


_ 1 «le(yj + f's) - «le<-'> 

c(-l) 


Y, 

1 2 c(-l) 


,2 (A2 + B2^ 


kk“ = 


4A2B2 


*2^= 


X 


1 SleW * b ) - «le<*> 


2 2 


Y _^2L:Jl 

2 2 c(l) 


c(l) 


(B17) 




If c(-l) = c(l) = 0 (i.e., the taper ratio is zero), the limiting values for kj and k^ 

must be used: 
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lim = 

c(-l)-0 ^ 


lim 

c(-l)-0 


- <=(yv) t's - 


T 2 

where d(y ) is the x-location of the midchord at y . Similarly, lim kg =1. It 
\ V c(+l)-*0 

follows that 


k'l In k\ - k'^ + kg In k^ - k^ - -2 


(B18) 


for a symmetrical wing in which c(l) — 0 and c(-l) — 0. 

The spanwise integration of the second and third terms in equation (B15) can be 
accomplished by Multhopp's method of integration. Substitution of equation (B12) for Lj 
results in 


-2n w 


2 " -h I ll ^3^ 


v)‘‘ 


N-1 . 

'--If 

2' J^o 


qjh) I, 


slgn(y^-q) 


In kM 


dr? 




Multhopp' s method of integration can be applied by using the trigonometric interpolation 
formula of the type 


m-1 


m 


V”> = ^ 21 (Vj)„ I ='* 

m-1 ^=1 

n= — — - 


sin X6 
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where 6 = and m is the number of spanwise integration stations over the 

n 2 m + 1 

whole wing. Thus, 


N-1 


cij(i?n) 




G, 


,dk’ 


sign (y^ - n„) g qj(yJ >” >''nn 


(B19) 


where the prime on the summation sign means that the term with n = is omitted from 
the summation and 


dkV k' 


Ml Ml 


dV 2A^B^ 




X(A - B)^ + - AB 


dx , [; . 2^.^ (A _ b,2] s 

dV 8k'.„A2B2 ^ J dr; 


. 1 ^le(yir) + H^v) ^sl - hei^n] 


''N 


{%) 


2 c 


A^ = (1 - X)^ + 


b2 = x2+y2 


dX 


- ^le(^n) ^ (jleN ~ «^s) ' ^e("n) ' 


> (B20) 


dr; 2 c(n^'^ 


2 c‘ 


’(^n) 


dr; 2 '|c(r;„) 


2 _ (A + B)^ - 1 
W1 4AB 


1 (y,i - \) ^e(^n) - 

V ^ 


® '"W 


G. 


cos j<^ + cos (j + 1) <p 
j('^s) = sin <t>g 


cos 0 + j sin + (j + 1) sin (j + 1) 


J 
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Also 


b 


vv 


m + 1 
4 sin 


b 


vn 


°‘”^n 

(m + 1) ^cos - cos 


= 0 




^|n - i/| odd 
^|n - v\ even 




(B21) 


hpu = ^s) 

_ ^^s ^sl 

jTl J 




Uil) 


(B22) 


Combination of equations (B16) and (B19) gives the total induced downwash: 


N-l C 

-2’'w(0g,y^)= J Jq.(y^) 


3=0^ 




V* o V* 

-z -sf ‘""’J-Z "m.’jW 


3=0 ' ' 


(y^) 


v' sin e„ G. f(p\ dk' 

y /3 £_li_^_^lnk' 

m + 1 drj ^ 




Jim/ 


(B23) 
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IMPROVEMENTS IN LANGLEY PROGRAM A0313 


The improvement in the logarithmic -singularity correction term presented in this 
report has been implemented into the analysis version of the computer program described 
in the supplement to reference 7 (Langley program A0313). In addition, 

(1) The program has been restructured into an overlay arrangement (core require- 
ments of 5IOOO3 words on the Control Data 6600 computer system). 

(2) The chordal loading functions of reference 5 are employed, 

(3) The test to set the values for N and m is utilized before a solution is begun. 

(4) Damping in roll, damping in pitch, and side -edge suction -force computations 
have been added. In order to access these items, the format of the third data input card 
has been changed to 5F6.0, F6.2, and 4F6.0 with the last three fields set aside to receive 
the codes needed to commence these computations. The codes are as follows: 


(a) In columns 43 to 48, a 1 causes the damping-in-roll stability derivative 
C, to be computed; a 0 indicates that C, is not desired. 

tp tp 


^mq 

that 


(b) In columns 49 to 54, a 1 causes the damping -in -pitch stability derivative 
and the lift coefficient due to pitch rate Cj to be computed; a 0 indicates 

q 

Cjn and Cj^ are not required. 


(c) In columns 55 to 60, a 1 causes the side -edge suction force to be computed 
and the aerodynamic characteristics as a function of angle of attack to be listed; 
a 0 indicates that neither the side -edge suction force nor a listing of the aerodynamic 
characteristics is required. 


Note that roll-rate and pitch-rate stability derivatives cannot be computed simultaneously. 
Also, the side-edge suction force should not be computed concurrently with roll-rate or 
pitch-rate derivatives. 

This version of the computer program is available from COSMIC, 112 Barrow Hall, 
University of Georgia, Athens, GA 30602. 
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